Global properties of locally spatially homogeneous 
cosmological models with matter 



Alan D. Rendall* 

Max-Planck-Institut fur Astrophysik 
Karl-Schwarzschild-Str. 1 
Postfach 1523 
85740 Garching 
Germany 



Abstract 

The existence and nature of singularities in locally spatially homogeneous solutions of 
the Einstein equations coupled to various phenomenological matter models is investigated. 
It is shown that, under certain reasonable assumptions on the matter, there are no singu- 
larities in an expanding phase of the evolution and that unless the spacetime is empty a 
contracting phase always ends in a singularity where at least one scalar invariant of the 
curvature diverges uniformly. The class of matter models treated includes perfect fluids, 
mixtures of non-interacting perfect fluids and collisionless matter. 



1. Introduction 

When trying to understand the global properties of solutions of the Einstein equa- 
tions, a tractable starting point is the study of solutions with high symmetry. Assuming 
that there is a symmetry group with three-dimensional spacelike orbits leads to the class 
of spatially homogeneous spacetimes. More generally, it is possible to consider space- 
times which are locally isometric to spatially homogeneous ones. For these spacetimes the 
Einstein equations reduce to ordinary differential equations, an enormous simplification. 
At the same time these spacetimes are general enough to display a variety of interesting 
dynamical behaviour. 
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The purpose of this paper is to study the global properties of locally spatially ho- 
mogeneous cosmological models which admit a compact Cauchy surface. Attention will 
be confined to certain phenomenological matter models. No attempt will be made to de- 
termine the details of the evolution of these spacetimes; the aim is rather to answer two 
fundamental questions. The first is whether singularities can occur in an expanding phase 
of the evolution of the spacetime and the second is whether a singularity which occurs in a 
contracting phase must be a curvature singularity. A review of results on these questions 
in the case where the matter is a perfect fluid can be found in [6] . 

The class of cosmological models to be considered will now be defined. They all admit 
Cauchy surfaces and have a matter content which is such that the Cauchy problem for 
the Einstein equations coupled to the equations describing the matter model is well-posed. 
Hence it is enough to define the class of initial data which gives rise to them. The initial 
data consists of a Riemannian metric g a b, a symmetric tensor k a b and some matter data, 
denoted collectively by Fq, on a three-dimensional manifold M. 

Definition An initial data set (g a b, kab, F ) for the Einstein-matter equations is called 
locally homogeneous if the naturally associated data set on the universal covering manifold 
M is homogeneous i.e. invariant under a transitive group action. 

The terminology is motivated by the following fact. A Riemannian manifold (M, g) is 
called locally homogeneous if given any two points x,y G M there exists an isometry <p 
of a neighbourhood of x onto a neighbourhood of y with <p{x) = y. Singer's theorem 
[14] says that a complete Riemannian manifold M is locally homogeneous if and only if 
its universal cover is homogeneous. It is left open here whether an analogue of Singer's 
theorem holds for initial data sets for the Einstein equations. The spacetimes considered 
in the following will be Cauchy developments of locally homogeneous initial data sets on 
compact manifolds. It is convenient for many purposes to work on the universal covering 
space M. In the notation no distinction will be made between objects on M and their 
pull-backs to M. 

The universal cover of the given spacetime may not be spatially homogeneous but it 
can be extended to be spatially homogeneous. (This follows from the fact that the equations 
describing the evolution of the geometry and the matter fields on the universal cover can be 
written in a form where there is no explicit dependence on the spatial variables.) It will be 
assumed that such an extension has been made. Then the universal cover has a preferred 
foliation by orbits. Each leaf of this foliation has constant mean curvature. There is an 
induced foliation of the original spacetime by constant mean curvature hyper surfaces. It 
is then topologically of the form M x /, where I is an interval. A time coordinate can be 
defined by requiring that it be zero on the initial hypersurface and that it be an arc length 
parameter increasing towards the future on the timelike geodesies normal to the initial 
hypersurface. This time coordinate is constant on each leaf of the preferred foliation. A 
time coordinate of this kind will be called Gaussian. Let the mean curvature of the leaf 
with time coordinate t be denoted by H{t). In the following the time orientation will be 
chosen so that the mean curvature of the initial hypersurface t = is non-positive. In 
other words the model is initially non-contracting. 

The spatially homogeneous spacetimes can be divided into two classes, the Bianchi 
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models and the Kantowski-Sachs models. (For general information on this see [16].) Most 
of this paper is concerned with the Bianchi case. The Kantowski-Sachs case is discussed 
in Section 5. For Bianchi models the universal cover M can be identified with a Lie group 
G. Tensors on G will be described in terms of frame components in a left invariant frame. 
The components of the induced metric and second fundamental form of the homogeneous 
hyper surfaces will be denoted by gij and kij respectively. They are functions of t. The 
mean curvature is given by H = g^kij. 

The matter models to be considered will be defined in terms of some general properties. 
As usual T af3 denotes the energy-momentum tensor. When a specific matter model has 
been chosen T af3 will be a functional of some matter variables, denoted collectively by F, 
and the spacetime metric g a p. In the following another quantity N a (called the particle 
current density) will be required. It is also assumed to be a functional of F and g a p. Now 
various properties which will be assumed at appropriate points will be listed. 

(1) T af3 V a Wp > for all future-pointing timelike vectors V a and W a (dominant energy 
condition) 

(2) T a/3 (g a p + V a Vp) > for all unit timelike vectors V a (non-negative sum pressures 
condition) 

(3) for any F and g a p the conditions V a iV a = and V a T a/3 = are satisfied (conservation 
conditions) 

(4) for any F and g a p the vector N a is future-pointing timelike or zero 

(5) for any constant C\ > there exists a positive constant C2 such that for any F and 
g a p with —N a N a < C\ and any timelike vector V a the following inequality holds: 

T a/3 V a Vp > C 2 {N a V a ) 2 

(6) for any constant C\ > there exists a positive constant C2 < 1 such that for any F 
and g a p with —N a N a < C\ and any unit timelike vector V a 

(g a p + V a Vp)T a ? < 3C 2 T^V a Vp 

(7) if a solution with the given symmetry of the Einstein equations coupled to the given 
matter model is such that the time coordinate defined above takes all values in the interval 
(ti, t 2 ) , if it is not possible to extend the spacetime so as to make this interval longer and 
if t\ or t 2 is finite then H(t) is unbounded in a neighbourhood of t\ or t 2 respectively. 

(8) for any constant C\ > there exists a constant C 2 > such that T a pT a/3 < C\ implies 
-N a N« < C 2 

Some comments will now be made concerning the physical motivation of some of these 
conditions. If a given type of matter can be considered as being made up of particles then 



3 



a particle current density iV a is defined. If the particles have positive rest mass then this 
vector is future pointing timelike or zero as required by condition (4). If the particles are 
massless then this condition is still satisfied except for very special types of matter where 
N a might be null. If particles cannot be created or destroyed then iV a is divergenceless as 
required in condition (3). It is not easy to give an intuitive interpretation of conditions (5) 
and (6). The meaning of (5) is roughly as follows. If matter is observed from a boosted 
frame then the particle density is multiplied by a 7-factor arising from the effect of Lorentz 
contraction on the volume element. The observed energy density is also affected in this way 
but picks up an additional 7-factor. Hence when a given matter distribution is considered 
from a boosted frame the multiplicative factor in the observed energy density behaves like 
the square of that in the observed particle density. As for condition (6) , it will be seen in 
Section 4 that for a perfect fluid it is related to the condition that the speed of sound should 
be bounded away from the speed of light. The given symmetry referred to in condition (7) 
will be Bianchi symmetry or Kantowski-Sachs symmetry, according to the context. 

In the following the two fundamental questions formulated above will be answered for 
locally spatially homogeneous spacetimes which are spatially compact and where the mat- 
ter content satisfies (some of) the conditions (l)-(8). In fact it will also be necessary when 
examining the issue of curvature singularities to assume that the spacetimes considered 
are not empty since in the vacuum case it is much more difficult to separate spacetimes 
with curvature singularities from spacetimes with Cauchy horizons. The question of which 
spatially compact, locally spatially homogeneous vacuum spacetimes have Cauchy horizons 
has been investigated in [3]. In the following charged matter and field-theoretic matter 
models (e.g. Yang-Mills, cr-models) are not considered. It would be of interest to know 
whether similar results hold in those cases. 

The paper is organized as follows. Section 2 contains results on geodesic completeness 
in an expanding phase while Section 3 is concerned with curvature singularities in a con- 
tracting phase. In both cases the hypotheses on the matter model are of a general nature. 
Section 4 is concerned with verifying that these properties hold for various specific matter 
models, namely perfect fluids, non-interacting mixtures of perfect fluids and collisionless 
matter. As already mentioned, the last section is devoted to Kantowski-Sachs models. 



2. Geodesic completeness 

In this section it will be shown that under general conditions there can be no singu- 
larity during an expanding phase of a locally spatially homogeneous spacetime. For this 
statement the existence of a compact Cauchy surface is irrelevant. An expanding phase 
means by definition a time interval (£1,^2) where H < 0. The statement that there are 
no singularities means that either the solution can be extended to an interval (£1, t' 2 ) with 
t' 2 > £2 or that it is future complete (i.e. that all inextendible causal geodesies are complete 
in the future direction). This statement is a consequence of the following theorem. 
Theorem 2.1 Let a spatially locally homogeneous solution of the Einstein equations 
coupled to a matter model satisfying the conditions (1), (2) and (7) be given. Suppose 
that a Gaussian time coordinate defined on this spacetime takes all values in the interval 
(£1,^2) and that the spacetime cannot be extended so as to enlarge this interval. Suppose 
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further that the mean curvature H satisfies limsup^ t2 H(t) < oo. Then t 2 = oo and the 
spacetime is future complete. 

Proof The proof given here is for Bianchi models; for Kantowski-Sachs models see Section 
5. Since spatial compactness plays no role we may assume without loss of generality that 
the spacetime is simply connected. Now some of the Einstein equations for Bianchi models 
are needed, namely the Hamiltonian constraint 

R - kijk ij + H 2 = 16nT 00 (2.1) 

and the evolution equation 

d t H = R + H 2 - 12nT 00 + Ang l3 T lJ . (2.2) 

Here R denotes the scalar curvature of the spatial metric. Combining these two equations 
gives 

d t H = kijk ij + MT 00 + 9ij T^). (2.3) 

Now kijk^ > \H 2 and so equation (2.3) together with conditions (1) and (2) implies 
that dtH > \H 2 . In particular this shows that H is non-decreasing. Thus, under the 
hypotheses of the theorem, the fact that t 2 = oo follows immediately from condition 
(7). The inequality also shows that if H is ever positive it must become infinite in finite 
time. Since it is already known that t 2 = 00 it follows that H < 0. There are now two 
possibilities: either H vanishes at some time or it is everywhere negative. 

Suppose that H = at some time to- Since H is everywhere non-positive it can be 
concluded that dtH is also zero there. From equation (2.3) it follows that kij = and 
T 00 = for t = to- The dominant energy condition implies that T a @ = at t = to an d 
hence (see [10], p. 94) everywhere. Thus in this case the solution must be a vacuum 
solution. Under these conditions the Einstein evolution equations imply that the Ricci 
tensor i?^ of is zero at to- Hence the metric gij (to) is flat. Uniqueness in the Cauchy 
problem for the vacuum Einstein equations now shows that in the first case the spacetime 
is obtained by an identification of (the whole of) Minkowski space and future completeness 
is clear. 

Now consider the second case. There the universe expands for ever. It has been shown 
by Lin and Wald [11] that for Bianchi IX models conditions (1) and (2) imply that the 
universe cannot expand for ever. Thus under the hypotheses of the theorem the Bianchi 
type cannot be IX. It can then be concluded that R < 0. (For a simple proof of this see 
[15].) Using this in (2.1) gives 

hjk ij < H 2 (2.4) 
The eigenvalues of k^ with respect to g^ are the solutions of 

det(fc y - \ gi j) = (2.5) 

Call them Ai, A2 and A3. Then 

Ai + A 2 + A 3 = H, (2.6) 
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and 

\j + \ 2 2 + \ 2 S = kijk^ . (2.7) 

Define pi = Xi /H so that 

Pi+P2+P3<l- 

Let K denote the subset of R 3 where the relations (2.8) are satisfied. The minimum value 
of pi on K will now be determined. By symmetry this will also be the minimum value of 
P2 and ps- The minimum must occur on the boundary of K and the boundary is given by 
replacing the inequality in (2.8) by equality. Since K is convex the minimum must occur 
when two of the eigenvalues are equal and using this its value can be computed to be —1/3. 
It follows that for any vector x l tangent to the hypersurfaces of homogeneity the following 
inequality holds 

k ij x i x j < -JftyyxV. (2.9) 

Consider now a future directed causal geodesic and let q l be the frame components of 
the projection of its tangent vector onto the hypersurfaces of constant t. Then 



Using the inequality d t H > ±H once more shows that 

H(t) > -3/(C + t) 
for some positive constant C. Combining (2.10) and (2.11) gives 



Hence 



< -iHgijjq? 



j t (\ogg ij q i q i)<2/(C + t). 

(gijqWr^ycii + t)- 1 . 



(2.10) 

(2.11) 

(2.12) 
(2.13) 



The affine parameter length of the geodesic up to a certain point is given by the integral 
of (e + gijq l q-')~ 1 / 2 with respect to t, where e is 1 for a timelike geodesic parametrized 
by arc length and zero for a null geodesic. Future geodesic completeness is equivalent to 
the divergence of this integral as the upper limit tends to infinity. This is guaranteed by 
(2.13). 



3. Curvature singularities 

The subject of this section is a contracting phase of the universe i.e. a time interval 
(ti,t2) where H > 0. Only Bianchi models are considered here - analogous results for 
spacetimes with Kantowski- Sachs symmetry are presented in Section 5. It was shown in 
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the last section that conditions (1) and (2) imply that H is non-decreasing. It follows that 
any extension of the solution to an interval (ti, t' 2 ) with t' 2 > ti must also be a contracting 
phase in this sense. Extend the time interval as far as possible to the future and call the 
result (ti,£*). It was also remarked in the last section that if H is positive somewhere 
it must become infinite in finite time. Hence t* < oo. In this section, in contrast to the 
previous section, the spatial compactness of the model plays an important role. Also all 
the conditions (l)-(7) will be used. Writing out the equation V a iV a = explicitly for a 
Bianchi model gives 

d t ((det g^N ) = -C^N^detg) 1 / 2 (3.1) 

Here Cj fc are the structure constants of the Lie algebra corresponding to the Lie group 
being considered. This equation should be thought of as being defined on the universal 
covering spacetime where an invariant frame exists. In the Bianchi types of class A where, 
by definition, C|- = the right hand side of (3.1) vanishes identically and the quantity 
(det g) 1 / 2 ^ is conserved. It will now be shown that the assumption of spatial compactness 
is enough to guarantee the vanishing of the right hand side of (3.1), even for class B models. 
To see this let X be any vector on the Lie group G which is left invariant and which projects 
to a smooth vector field X on M. The divergence of X is constant because it is invariant 
and the divergence of X is the same constant. However the integral of div X over M must 
be zero by Stokes' theorem. It follows that X has divergence zero. In terms of components 
this divergence is just C\aXK Now is at each time a vector field whose pull-back to 
the universal covering manifold has the properties of the vector field X just discussed. It 
follows that the right hand side of (3.1) is zero as claimed. 

The aim now is to show that if conditions (l)-(7) hold and if N a does not vanish then 
the singularity which occurs as t — > £* is a curvature singularity. Because of the conserved 
quantity which has been found it suffices to assume that N a does not vanish on the initial 
hypersurface to ensure that it does not vanish anywhere. 

Theorem 3.1 Let a solution with Bianchi symmetry of the Einstein equations coupled to a 
matter model satisfying conditions (l)-(7) and admitting a compact Cauchy hypersurface 
be given. Suppose that a Gaussian time coordinate defined on this spacetime takes all 
values in the interval (£i,£*) and that the spacetime cannot be extended so as to enlarge 
this interval. Suppose further that H is positive somewhere and that the particle current 
density A^ Q does not vanish. Then the mass density measured by any observer is unbounded 
as t — > £*. If in addition condition (8) is satisfied then limsup^^ G af3 G a p(t) = oo, where 
G a p is the Einstein tensor of g a p. 

Proof The first thing to be done is to obtain some control on the rate at which H(t) blows 
up as t — > £*. (That it does blow up as t — > £* follows from condition (7).) It was already 
mentioned in the previous section that dtH > \H 2 . Integrating this inequality gives 

H(t) < (3.2) 

On the other hand (2.2) and the dominant energy condition imply that 

d t H <R + H 2 (3.3) 
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Let R+(t) = max(R(t), 0). Then d t H < R+ + H 2 . For t close to £* we have H > 1. Thus 
d t H < R + H + H 2 and integrating this inequality gives the estimate 

H(t) > (£* - t)' 1 exp (- R + (t')dt'^j (3.4) 

Here the right hand side is defined to be zero if the integral is infinite. If the Bianchi type 
is not IX then R + = and (3.4) simplifies to 

H^yiU-t)- 1 (3.5) 

From these estimates for H it is possible to get estimates for det g since d/dt (log det g) = 
-2H. It follows from (3.2) that 

detg > d{U -tf (3.6) 

for some positive constant C\. If J** R + (t)dt < oo for any t < then (3.4) implies that 

det g(t) < C 2 (U - t) 2 (3.7) 

The integrability condition on the scalar curvature is obviously satisfied when the Bianchi 
type is not IX. In type IX a more detailed analysis is necessary to see whether it holds. 
When (3.7) holds the fact that iV°(det g) 1 ^ 2 is conserved shows that 

iV° > C 3 (U -t)' 1 (3.8) 

Applying condition (5) to the unit normal vector to the homogeneous hypersurfaces gives 

T 00 > CN 2 (3.9) 

as long as N a N a is bounded. 

Suppose now that the Bianchi type is not IX. It will be shown that the assumption 
that N a N a is bounded leads to a contradiction in this case. If this quantity is bounded 
then (3.8) and (3.9) imply that there exists a constant C > such that 

T 00 >C(U-t)~ 2 (3.10) 

Then the Hamiltonian constraint (2.1) and the fact that R < when the Bianchi type is 
not IX give an estimate of the form 

hjk lJ < H 2 (l-7]), v >0 (3.11) 

Taking a linear combination of (2.1) and (2.2) gives 

d t H = a(k i:j k ij - H 2 ) + (1 - a)R + H 2 - 4tt(3 - 4a)T 00 + 4n gij T ij (3.12) 



Because of (6) it is possible to choose a > so that the sum of the last two terms in this 
inequality is not positive. Hence 

d t H < (1 - ar])H 2 (3.13) 
This allows the inequality (3.8) to be improved to 

N° >C(t*-t)T=^ (3.14) 

From (3.9) it follows that 

T 00 > C{U -t)T=*v (3.15) 

This contradicts the Hamiltonian constraint. The conclusion is that if the Bianchi type is 
not IX then N a N a must be unbounded. Consider now any observer, represented mathe- 
matically by a timelike curve with unit tangent vector t a . Then the mass density measured 
by this observer is 

-N a t a > ^f-N a N« (3.16) 

This gives the first conclusion of the theorem in the case that the Bianchi type is not IX. 
The other conclusion follows immediately from the fact that N a N a is unbounded. 

Consider now the Bianchi IX case. It is a standard fact (see e.g. [15]) that the 
dimensionless quantity (det g) x ^R is bounded above. Hence R < C(det g)~ x ^. Suppose 
for a moment that 

detg >C(U -t) (3.17) 

for some positive constant C. Then R + < C(t* — t) -1 / 3 . This implies the integrability 
of R + and hence, by (3.4), that H(t) > C/(t* — t). It follows that det g tends to zero as 
t — > £*. If, on the other hand, the inequality (3.17) does not hold for any constant C then 
there must exist a sequence of times t n with t n — > £* such that 

det<7(£ n ) < C(U -t n ) (3.18) 

Hence, in particular, &etg(t n ) — > 0. Since det g is monotonically decreasing near the 
singularity it follows that det g — > in this case too. 

It will now be shown that the assumption that N a N a is bounded leads to a contra- 
diction. This suffices to complete the proof of the theorem since from that point on the 
argument is identical to that given above for the case where the Bianchi type is not IX. In 
the equation (2.2) the combination R — 127rT 00 + Ang^T 1 ! occurs. If N a N a is bounded 
then condition (6) allows the sum of the second and third terms to be bounded from above 
by — CT 00 for some positive constant C. Also condition (5) gives 

T 00 > C(iV ) 2 = C"(det g)- 1 (3.19) 

Thus as det g tends to zero CT 00 becomes larger than R and the above combination of 
terms from (2.2) is negative. Hence d t H < H 2 . It follows that the inequalities (3.7), (3.8) 
and (3.10) remain true in the Bianchi IX case. The inequality (3.11) can also be obtained if 
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we replace the fact that R < used before by the fact that sufficiently near the singularity 
R < CT 00 for any given positive constant C. This can also be used to control the term 
involving R in (3.12) using a small constant multiplied by T 00 . This is enough to prove 
(3.13)-(3.15) in the Bianchi IX case. The same trick of dominating R by a small amount 
of T 00 shows that (3.15) contradicts the Hamiltonian constraint in this case also and this 
completes the proof of the theorem. 

Remarks 1. When applied to a general perfect fluid solution of Bianchi type IX this 
theorem is a non-trivial generalization of the results given in [6] . 

2. In [12] it was claimed that a singularity in a contracting phase of a Bianchi class A non- 
vacuum solution of the Vlasov-Einstein system must be a curvature singularity. The claim 
is not justified by the arguments given there but it is true. It follows from Theorem 3.1 
above and the fact, proved in the next section, that the Vlasov equation satisfies conditions 
(l)-(8). 

3. For solutions where iV Q is orthogonal to the hypersurfaces t =const. ^/—N a N OL = N° 
and so in that case the theorem remains true if conditions (5) and (6) are dropped from 
the assumptions. 



4. Matter models 

Before specific matter models are studied some general remarks will be made concern- 
ing condition (7) in spacetimes with Bianchi symmetry. Suppose that a solution is given 
on an interval (tijtz), that t 2 is finite and that H is bounded on a neighbourhood of t 2 . 
Then det g and its inverse are bounded on a neighbourhood of t 2 . It is a standard fact for 
Bianchi models that under these conditions the scalar curvature R is bounded above (see 
e. g. [15]). Next the Hamiltonian constraint shows that k^kij is bounded. An argument 
given in [12] (p. 88-89) then shows that and k^ are bounded in a neighbourhood of t 2 . 
Analogous statements hold in the case that the roles of ti and t 2 are interchanged. What 
is needed to check condition (7) is that when a solution of the Einstein-matter equations 
is given on some open interval and when gij, kij and (det g)~ x are bounded near an end- 
point of this interval, then the solution extends through this endpoint to a strictly longer 
interval. Call this condition (7'). 

It is now time to look at some examples. The matter model which has been studied 
most extensively in the context of Bianchi models is the perfect fluid. It will now be inves- 
tigated under which conditions on the equation of state a perfect fluid satisfies conditions 
(l)-(8). The equation of state is a relation p = f(p) between energy density and pressure. 
The energy-momentum tensor is T a ^ = pu a u^ + (p + p)g a @ The equation of state will be 
assumed to satisfy the following conditions 

(i) / is a continuous function from [0, 00) to the non-negative real numbers with /(0) = 
which is C 1 for p > 

(ii) 0<df/dp<l for each p > 

Condition (1) is satisfied if \p\ < p ([10], p. 91) while condition (2) is equivalent to 
p > 0. Hence both of these conditions are consequences of (i) and (ii). The quantity A^ a 
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is denned to be ru a , where 



p 



exp / (p' + ZtpOrV (4.1) 



This definition ensures that the Euler equations V Q T a/3 = imply V a -?V a = and hence 
that (3) is satisfied. Condition (4) obviously holds. A sufficient condition for (5) to be true 
is that an upper bound on r implies an inequality of the form p > Cr 2 . In fact this follows 
from (4.1) and the dominant energy condition. Condition (6) represents a restriction on 
the equation of state which goes beyond (i) and (ii), namely that when p is less than some 
given constant C\ there is a constant < 1 such that p < C^p- It will be expressed by 
saying that the fluid is not asymptotically stiff at low densities. Condition (8) obviously 
holds. 

The one condition which remains to be verified for a perfect fluid is (7 ; ). In order to 
do this it is useful to write the Euler equation V a T a ^ = in terms of the density p and 
the spatial components v l of the four- velocity (expressed in an invariant frame). Then no 
spatial derivatives occur and a system of ordinary differential equations is obtained. If we 
denote the matter variables p and u l collectively by / then the equations take the form 

A(f,g ij )df/dt = F(f,g ij ,k ij ), (4.2) 

where the matrix A and the vector F depend smoothly on their arguments. As long as 
p > the determinant of A is non- vanishing. Hence the evolution equation can be written 
as df/dt = A~ X F. The right hand side of this latter equation is regular so long as detg 
and p are bounded away from zero. This equation can be combined with the evolution 
equations for the geometry, which in the homogeneous case are 

d t 9ij = -Zkij (4.3) 
d t k i:i = Rtj + {g lm k lm )k %J - 2k u k) - 8nT tJ - 4nT 00 g tJ + 4n(g lm T lm )g i:j (4.4) 

The Ricci tensor Kij is a rational function of g^ whose exact form depends on the Bianchi 
type. Taking the equation for df/dt together with (4.3) and (4.4) gives a system of ordinary 
differential equations for /, g^ and k^. The coefficients in these equations are regular as 
long as detg and p are strictly positive. Standard results on existence and uniqueness 
for ordinary differential equations (see e.g. [9]) show that given initial data (p, u\ g^, k^) 
with p > there exists a unique solution of these equations on some time interval with the 
given initial values. Moreover this solution can be extended as long as u l , gij, kij, p, p~ x 
and (detg) -1 remain bounded. 

To verify condition (7') it remains to show that if a solution of the Einstein-Euler 
system is given on some time interval and if g^, k^ and (detg) -1 are bounded then the 
relevant matter quantities are bounded. The conservation law for shows that ru° and 
its inverse are bounded. Hence it can be concluded that r is bounded above. For a perfect 
fluid condition (6) implies that for r < C\ an inequality of the form p > C2r 2-e holds, for 
some e > 0. Now 

T 00 > p(u ) 2 > Cr 2 -'(u° f 
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The Hamiltonian constraint shows that T 00 is bounded and so it follows that r _1 is 
bounded. Combining this with the boundedness of ru° shows that u is bounded. This 
completes the verification of condition (7') for a perfect fluid. 

The above information concerning perfect fluids will now be collected. 
Proposition 4.1 Consider a perfect fluid with equation of state p = f(p) satisfying 
conditions (i) and (ii). Then conditions (l)-(5)and (8) are satisfied. If, in addition, the 
fluid is not asymptotically stiff at low densities then conditions (6) and (7) are also satisfied. 

It follows from this proposition that Theorem 3.1 applies to any perfect fluid which satisfies 
the conditions (i) and (ii) and which is not asymptotically stiff at low densities. 

Another kind of matter model which has been popular in cosmology is that of n 
(usually two) non-interacting perfect fluids [4]. The energy-momentum tensor in this case 
is of the form 

where each summand has the form of the energy-momentum tensor of a perfect fluid 
with its own four- velocity u%y density and pressure p^y The latter are related by 
equations of state pu\ = fu)(p(i))- The condition that the fluids are non-interacting is 
expressed by the equations V a T^ = for each i. If we define by putting an index 
i on all the quantities in (4.1) and let = r^ufy then it follows that V a iV^ = for 
each i. In particular the sum N a of the Nfy has zero divergence. Before showing under 
which assumptions a mixture of perfect fluids satisfies conditions (l)-(8), a more general 
result will be given on matter models which are composites of two others. 
Lemma 4.1 Let two matter models be given with matter variables and F( 2 ), particle 
current densities iV^ and iV^ and energy-momentum tensors and T^. Define a 
new (composite) matter model by taking as matter variables F = (i^i), ^(2)) an d defining 
N a and T a(i as the sums of the corresponding quantities for the individual matter models. 
If the original matter models satisfy conditions (l)-(6) then the composite model does so 
too. 

Proof It is obvious that the composite matter model satisfies conditions (l)-(4). Next 
note that -N a N a > -NfaN^ - N{* 2) n£\ so that the boundedness of N a N a implies 

that of NfaNa for i = 1,2. It follows that condition (6) is satisfied and condition (5) is 
a consequence of the following computation: 

TapV a VP = T^V a V p + T { ^V a V p 

> C 2 [{N^V a ) 2 + (N^V a ) 2 ] (4-7) 

> \C 2 {N«V a ) 2 



Consider now a non-interacting mixture of perfect fluids and suppose that each individual 
fluid satisfies the hypotheses of Proposition 4.1 including the condition that it should not be 
asymptotically stiff. It then follows from Lemma 4.1 that the mixture satisfies conditions 
(l)-(6). It will now be investigated whether conditions (7) and (8) also hold. The evolution 
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equations for a mixture of non-interacting fluids can be written in a similar way to those 
for a single fluid. Thus to verify condition (7) in this case it suffices to check that the fluid 
quantities remain finite in a given regular spacetime. However each fluid satisfies the Euler 
equations separately and so the argument given above for a single fluid applies without 
change. (It must, of course, be assumed for this that none of the fluids are asymptotically 
stiff at low densities.) A straightforward calculation shows that 

T a(5 T a0 >(£p {i) ? (4-8) 

(0 

Hence if T a/3 T a/ 3 is bounded all the p^ are bounded. Moreover, the boundedness of 



T af3 T a p implies that the quantities p^p^(u?^Ua) 2 are bounded for each It follows 



that N^N^ is bounded for each Thus condition (8) is also satisfied and the following 
analogue of Proposition 4.1 is obtained. 

Proposition 4.2 Consider a mixture of non-interacting perfect fluids, each of which sat- 
isfies the hypotheses of Proposition 4.1 including the condition that none of them be 
asymptotically stiff at low densities. Then conditions (l)-(8) are satisfied. 

The next matter model which will be considered is a collisionless gas described by the 
Vlasov equation. For an introduction to the Vlasov-Einstein system see [7] . Some general 
information on spatially homogeneous solutions can be found in [12]. In the following we 
consider a gas of particles of unit mass. The distribution function / is supposed to be 
a C 1 function with compact support in velocity space initially and these properties are 
preserved by the time evolution. The quantities N a and T af3 are defined in the usual way 
as certain weighted integrals of /. Conditions (l)-(4) are standard properties of solutions 
of the Vlasov-Einstein system. Condition (7) for this matter model is the statement of 
Lemma 2.2 of [12]. It remains to check conditions (5), (6) and (8) and to do this we may 
choose a frame whose timelike member is V a in order to do the calculation. Then the 
inequalities of (5) and (6) become f 00 > C 2 (N ) 2 and 5 tJ f ZJ < 3C 2 T 00 . (The hats here 
indicate the use of indices associated to an ortho normal frame.) 

= - f j f(p a )f(q a )p a q a /(p°q )dp 1 dp 2 dp 3 dq l dq 2 dq 3 
f f f(p a )f(q a )/(p\ G Wdp 2 dp*dq 1 dq 2 dq* 
f(p a )/p°dp 1 dp 2 dp s 



(4.9) 



> 



Hence 



N° = J f(p a )dp 1 dp 2 dp 3 

/ f(p a )p°dp 1 dp 2 dp 3 J ( / f(p a )/p°dp 1 dp 2 dp 3 



(4.10) 



< (T 00 ) 1 /2(_iV Q iV«)i/4 



13 



This shows that (5) holds. It follows directly from the definitions that the inequality of 
(6) holds with C 2 = 1/3 even without restricting N a N a to be bounded. Finally 



Proposition 4.3 Consider a collisionless gas of particles of unit mass for which the distri- 
bution function is C 1 and compactly supported in velocity space. Then conditions (l)-(8) 
are satisfied. 

Most of what was said in the last paragraph only used the fact that matter is described 
by a distribution function and the expression for the energy-momentum tensor in terms 
of this distribution function. Hence it has more general applicability to matter models 
based on kinetic theory. Consider the case of the Boltzmann equation. It follows from the 
above that conditions (l)-(6) and (8) are satisfied for the Boltzmann equation. The one 
condition which causes difficulty is (7). Recall furthermore that to verify (7) it is enough 
to verify (7'). For the Boltzmann equation this is a non-trivial task. Even in the case of 
spatially homogeneous solutions of the classical Boltzmann equation it is difficult to prove 
global existence and uniqueness theorems. (For details see [2] and references therein.) Just 
how difficult it is depends on the form of the collision kernel. For the special relativistic 
Boltzmann equation there appear to be no results in the literature on global existence of 
unique solutions for general homogeneous initial data. (There is a result for inhomogeneous 
data close to equilibrium [8].) 



5. Kantowski-Sachs spacetimes 

In this section a contracting phase of the universe will be considered as in Section 3 
but now the symmetry is of Kantowski-Sachs type rather than Bianchi type. It should 
be mentioned that for this symmetry type no singularities occur in an expanding phase 
if conditions (1), (2) and (7) hold. It can be shown as in the Bianchi case that if (7) 
holds either the solution can be extended to the future beyond the expanding phase or 
the solution exists globally in a Gaussian time coordinate and expands for ever. The 
impossibility of the second alternative under conditions (1) and (2) has been shown by 
Burnett [1]. 

Suppose now that a solution is defined on an interval (t±, £*) which has been maximally 
extended towards the future. Suppose further that this is a contracting phase. In a 
Kantowski-Sachs spacetime the metric can be written in the form 




(4.11) 



= -N a N- 



a 



ds 2 = -dt 2 + A 2 (t)dx 2 + B 2 (t)dE 2 



(5.1) 
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where dY? is the standard metric on the two-sphere. It is easy to check that the divergence 
of the second fundamental form of the hypersurfaces t=const. is zero. It follows from the 
momentum constraint that d/dt is an eigenvector of the energy- momentum tensor with 
eigenvalue T 00 . Hence 

T af3 T^ > (T 00 ) 2 (5.2) 

It can be shown that, as in the case of the Bianchi models, (det g) l / 2 N° is time independent. 
In this case det g denotes the determinant of the three-dimensional metric in the above 
coordinates, so that (det fir) 1 / 2 = AB 2 . The following theorem is an analogue of Theorem 
3.1 for spacetimes with Kantowski-Sachs symmetry. 

Theorem 5.1 Let a solution with Kantowski-Sachs symmetry of the Einstein equations 
coupled to a matter model satisfying conditions (l)-(7) and admitting a compact Cauchy 
hypersurface be given. Suppose that a Gaussian time coordinate defined on this spacetime 
takes all values in the interval (ii, £*) and that the spacetime cannot be extended so as to 
enlarge this interval. Suppose further that H is positive somewhere and that the particle 
current density iV a does not vanish. Then t = £* is a curvature singularity in the sense 
that some curvature invariant blows up uniformly as t — > £*. 

Proof Note first that, as a consequence of condition (7) H tends to infinity as t — > £*. A 
spacetime with Kantowski-Sachs symmetry is, in particular, spherically symmetric. It is 
thus possible to define the area radius r which, in terms of the coordinates used above, is 
just B. The gradient of r is everywhere timelike and its length squared is — B 2 . The mass 
function m is defined in a spherically symmetric spacetime by 

1 - 2m/r = V a rV a r (5.3) 

Thus in the Kantowski-Sachs case 

m=\B{l + B) 2 (5.4) 

In a spherically symmetric spacetime m is monotone increasing along the integral curves of 
W as long as the latter are timelike[l]. Hence m can be bounded from below in terms of a 
positive constant only depending on the initial data. When studying spherically symmetric 
spacetimes it is useful to consider quantities defined on the quotient of spacetime by the 
group action. This space inherits a 2-dimensional metric g a b with Gaussian curvature K 
and a 2-dimensional energy- momentum tensor T^. Let trT = g ab T a b. The Kretschmann 
scalar is given by 



R a ^ 5 R a ^ S = 4if 2 + 4r" 4 (V a rV a r) 2 + r" 4 (^V a rV a r - 2yrrtrT^ 
+ 167r 2 r- 2 (r ab - itrT(7 a6 ) (V 6 - ^trT^ 6 



(5.5) 



This formula also applies to spacetimes with plane symmetry and was used in [13] in that 
context. All terms except the last are manifestly non-negative. The last is also non- 
negative provided the energy-momentum tensor is diagonalizable. For Kantowski-Sachs 
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symmetry this is automatic. Hence the Kretschmann scalar can be bounded below by 
16m 2 /r 6 . Since it is known that m is bounded below by a positive constant it follows that 
if lim^^ r = the Kretschmann scalar will blow up as t — > t*. Thus to prove the theorem 
we can assume without loss of generality that r does not tend to zero at the singularity. 
The scalar curvature R of the t=const. hypersurfaces is simply B~ 2 . Thus under the 
present assumptions R is bounded above. Given this information (3.4) can be used as in 
the Bianchi case (when the type is not IX) to show that detg — > as t — > t*. By the 
conservation law N° tends to infinity in this limit. Condition (5), applied to the normal 
vector, shows that T 00 tends to infinity in this limit. In the Kantowski-Sachs case T 00 
is an eigenvalue of the energy-momentum tensor and so T a/3 T a( g — > oo as t — > t*. This 
completes the proof. 

It may be seen as unsatisfactory that this argument leaves open the nature of the 
curvature singularity in any given case. All that the above proof shows is that either 
R a ^ s R a p 1 s or G af3 G a p blows up as the singularity is approached. It will now be shown 
with the help of arguments of Collins [5] that the second of these quantities always blows 
up. According to the above proof, if B does not tend to zero as t — > £* then the statement 
is true. One of the field equations reads 

2BB + 1 + B 2 = -TlB 2 (5.6) 

If B tends to zero as t — > t* then B must be negative somewhere. Equation (5.6) then 
implies that it must stay negative. It follows from (5.6) that d/dt(B(l + B 2 )) > 0. In fact 
this is just the monotone property of the mass function mentioned earlier. It follows that, 
fixing some time to an d letting B = B(to): 

B (tf > -1 + B(*o)(l + B(t ) 2 )/B(t) (5.7) 

for any t > to. Since B — > it follows that B — > oo. Hence in particular for any e > 
the inequality e(B/B) 2 > 1/B 2 holds close to the singularity. The explicit form of the 
Hamiltonian constraint in this case is 

2AB B 2 1 00 

-Jb + ~B 2 + ~B 2 = MT (5 - 8) 

Hence 

This implies an inequality of the form 

H>C\B/B\ (5.10) 

In other words -d/dt (log det g) > -Cd/dt(logB) or d/dt(det g/B c ) < 0. It follows that 
det g tends to zero as t — > t*. Hence T 00 is unbounded and the same is true of G a ^G a f3. 

In order to be able to apply Theorem 5.1 to the matter models discussed in Section 
4, it is necessary to verify that condition (7) holds for these matter models in the case 
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of Kantowski-Sachs symmetry. Note first that (7) follows from condition (7') as in the 
Bianchi case. The proof is similar to that given for Bianchi symmetry at the beginning of 
Section 4. In the Kantowski-Sachs case it is not known a priori that the scalar curvature 
R is bounded above. However the evolution equation (2.3) for the mean curvature shows 
that on an interval where H is bounded the integral in time of kijk lJ is bounded. This 
is enough to apply the argument of [12] to obtain the boundedness of the metric. This, 
together with the boundedness of (det shows that R is in fact bounded and from that 
point on the argument proceeds as before. In the case of collisionless matter condition (7) 
was verified in [12]. For a perfect fluid of mixture of perfect fluids (7 ; ) can be checked as 
much as in the Bianchi case. The only difference is in the parametrization of the 4-velocity 
of each fluid. This necesarily lies in the (t, x) plane and so can be parametrized by its 
component in the x-direction. 
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sion of this paper. 
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